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Abstract 



Noether's theorem is reviewed with a particular focus on an intermediate step between global and local 
gauge and coordinate transformations, namely linear transformations. We rederive the well known result 
that global symmetry leads to charge conservation (Noether's first theorem), and show that linear symmetry 
allows for the current to be expressed as a four divergence. Local symmetry leads to identical conservation 
of the current and allows for the expression of the charge as two dimensional surface integral (Noether's 
\Q • second theorem). In the context of coordinate transformations, an additional step (Poincare symmetry) is 

| of physical interest and leads to the definition of the symmetric Belinfante stress-energy tensor, which is 

then shown to be identically zero in generally covariant first order theories. The intermediate step of linear 
symmetry turns out to be important in general relativity when the customary first order Lagrangian is used, 
which is covariant only under affine transformations. In addition, we derive explicitely the canonical stress- 
energy tensor in second order theories in its identically conserved form. Finally, we analyze the relations 
between the generators of local transformations, the corresponding currents and the Hamiltonian constraints. 



1 Introduction 

Our motivation originates in the following question: Why is it that in gravity theory, the stress-energy (pseudo)- 
tensors usually presented in literature are always in a form such that they are identically divergence free, while 
this is not the case for special relativistic theories? One could try to attribute this to the special role played by 
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the gravitational field, which is supposed to define the geometry of spacetime itself. But on the other hand, one 
can also take a less geometric point of view and forget about metrical structures and geometry in general. We 
are then left with a set of points x l (the spacetime manifold) and fields defined on this set, one of which being 
the gravitational field gn~. From this viewpoint, there is no particular reason why the stress-energy tensor for 
the gravitational field should have any essentially different properties from that of any other field. 

The answer to the above question is given by Noether's second theorem. While global symmetry of a 
Lagrangian theory leads to a conservation law (first Noether theorem), invariance of the theory under the same, 
but localized symmetry leads to identical conservation of the same current, in the sense that the charge can be 
written in the form of a surface integral. The exact formulation of the theorem can be found in the corresponding 
literature, in particular in Noether's original article pQ. In this paper, we will review what it means concretely 
in the case of gauge and in particular of coordinate transformations. The answer to our initial question is 
provided by the fact that gravitational theories are covariant under general coordinate transformations, while 
special relativistic theories are only Poincare invariant. 

Since the results presented in this paper were already known to Noether herself, there is no need to provide 
a large reference list. Our aim is, on one hand, to present the results in a mathematically simple form (as 
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opposed to fiber bundle descriptions) for the cases of physical interest, and, on the other hand, to highlight the 
relevance of the different degrees of locality of the symmetry. 

Our main focus lies on the stress-energy tensor and the conservation of energy and momentum, but in 
order to illustrate the procedure, and in particular the role of the intermediate step between global and local 
transformations, we start with a simple example of an internal gauge symmetry. In the first five sections of this 
paper, we assume that the Lagrangian C does not depend on second and higher derivatives of the fields. We 
refer to such theories as first order theories. Second order theories will be considered in section 6. 

The paper is structured as follows. In the remaining part of the Introduction, we derive Noether's theorem 
for first order non-abelian gauge theory. Then, in sections 2 to 4, we turn to coordinate transformations and the 
stress-energy tensor. We start with global translations and Poincarc symmetry in section 2, go over to affinc 
symmetry in section 3, and finally to general covariance in section 4, where, in each step, the consequences 
of the symmetry on the conservation law for the canonical as well as for the Bclinfante stress-energy tensor 
are discussed. In section 5, we apply those results to concrete theories, in particular to general relativity and 
to Einstein-Cartan theory. In section 6, we generalize the formalism to include theories based on Lagrangians 
containing second derivatives of the fields. Further, in section 7, we briefly analyze the relations between the 
Bclinfante and the Hilbert (or metric) stress-energy of the matter fields. Finally, in section 8, we analyze the 
relations between the generators of gauge and coordinate transformations, the corresponding currents and the 
first class Hamiltonian constraints of the theory. 

We begin by considering a Lagrangian C depending on a matter field ip and a gauge field Af , and assume 
that C is invariant under the following transformation 



F a 4- r a A 1 V 3 



Stp — —ie a a a ip. 



(1) 



The variation of C then reads 1 



dC SAt + ^f-SA 1 ^, 
dA? 4 fc '' 



Taking into account the field equations d£/dip = (<9£/<9V,i),i and dC/dAf, = (dC/dA^A^, we find 



= 



. dC , dC a 
- % W^ + JAt Ca " Al 



k,i 



e a + 



dC , dC a A „ . dC . 



dC 



Let us define the current by J l a 



cr a ip 



c a ^4j. In a first step, we assume invariance of C under 



dC 



global transformations, i.e., s" — 0. We then find the well known conservation law 

Ji 4 - o. 



(2) 



Next, we assume that the Lagrangian is invariant (in addition) under linear transformations, i.e., e a = s a (x), 



with e°j k = 0. This leads to the relation 



Ji _ ( 9C N 



(3) 



which expresses the current in the form of a divergence. In a last step, we assume local gauge invariance (i.e., 
general e Q (x)), which gives us a third relation 

dC 



OA 



= 0. 



(4) 



(»,fe) 



In order to shorten our expressions, we omit the terms in i/>, which arc similar to those in ip. 
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To summarize, global invariance leads to the conservation law for the current (and thus to charge conser- 
vation). Linear invariance allows us to write the current in the form Q, and finally, local symmetry leads to 
which tells us that the expression J3J is the divergence of an antisymmetric quantity, which as such has 
an identically vanishing divergence. Thus, in a locally invariant theory, the current can be written in a form 
such that it is identically conserved. In particular, this means that the charge (integration over a spacelike 
hypersurface) 

Q a = / J*jtei = J J°d 3 x (5) 

can be written in the form (greek indices from the middle of the alphabet, fj,,u . . . refer to the spatial components 
and run from 1 to 3) 

/ dC , / dC l2 

Qa = -JW, * k = ~fW» - d ^> (6) 

i.e., in the form of a two dimensional surface integral. Note that this would not be possible without the relation 
J2J. In the abelian case, equation (jfJJl is nothing but the well known expression of the charge in form of an 
integral of E over a closed surface. 

Similar results will be derived in the next sections for coordinate transformations. 



2 Translations and Poincare transformations 

We consider now a generic Lagrangian C = £(<p, <p j), where tp denotes collectively all the fields (including the 
gravitational field, whenever present). Under a global coordinate transformation x l — > x % — £ l , with constant 
all fields (scalar, vector, tensor, spinor) transform as 

&<P = Vtt, (7) 

where Sip denotes the change of (p at the point with the same coordinates, i.e., Sip(x) = <p'(x) — <p(x). If we 
assume invariance of the action under global coordinate translations, the Lagrangian must transform as scalar 
and we find 

SC = CiC = —Sip + - — Stp ti 
dip d<pi 



which leads to the conservation law 
with the canonical stress-energy tensor 



T l k ,i = 0, (8) 

t\ = £L<p tk - SIC. (9) 

Note that, although C will, in generally covariant theories, be a scalar density rather than a scalar, we will 
nevertheless stick to the definition JjJJ (thus including the factor \J—g into T % k ). 

Thus, our first result is that global translational symmetry leads to a conserved energy-stress tensor. It is 
not hard to see that every Lagrangian that does not depend explicitely on x % possesses that symmetry. 

Before we turn to linear transformations, we consider Poincare transformations, which are of particular 
physical interest. The fields now transform as 

Sip = ?<p,i + \e ik {Sv) ki , (10) 
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with £ l = a 1 + s l k x k (constant a 1 and e z k ) and e lk — e l l rj lk — —e kl . By (Sip) k i, we denote the action of the 
Lorentz group on the field in question. Explicitely, we have 

(Sip) hi = scalar (11) 

{SAi) ki = rjkiAi - ijuA k vector (12) 

(Shim) k i — Vkih im - rjuhkm + Vkm.hu - Vimhkl tensor, (13) 

and similar for contravariant or mixed tensors A % ,h lk ,h % k . Note that, since we are ultimately interested in 
generally covariant theories, we have to assume that spinors (spin 1/2) are described by fields transforming 
as scalars under Lorentz transformations 2 (since else, the generalization to general linear and diffeomorphism 
transformations would not be possible with finite dimensional representations). Such fields differ from true 
scalar fields (spin 0) by their behavior under local Lorentz gauge transformation, unrelated to the coordinate 
transformations we consider here. 

Assuming Poincare invariance of the action, and thus 8C = we find after some simple manipulations, 

apart from JHJ), the additional relation 

(r lk + \s mlk m )e lk = 0, 



where r = rj r l ; and 



S mik = -^-(SipY k . (14) 



Thus, since e is arbitrary (but antisymmetric), we get 

T [ik] + \ S mlk m = 0. (15) 

This result can be compared with J2J, but it is less strong since it determines only the form of the antisymmetric 
part or r lk . However, Ijl5(l can be used for a different purpose, namely the symmetrization of r lk . Indeed, we 
can define the so-called Belinfante tensor (see [5] and 0) 

T lk = r tk + I[5" ifcm - S mkl - S kmi ], mi (16) 

which obviously satisfies T %k ' i = 0, since the expression in brackets is antisymmetric in im. In other words, T lk 
differs from r lk only by a so-called relocalization term of the form C %mk m , with C vmk = — C" rllfc , i.e., a term 
whose divergence C tm m i vanishes identically and which leads only to two dimensional surface terms in the field 
momentum. Also, using (|15|) . it follows immediately that we have 

rpik rpki (17) 

Therefore, we can also define a spin current density in the form <j lk% = T tl x k — T lk x, satisfying <j lk% . = 0, which 
is not possible with the asymmetric tensor r lk . 

We will further discuss the Belinfante tensor later on. 



2 Spin 3/2 fields ipi carry an additional vector index and transform according to 1121 . 
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3 Linear transformations 



Let us also define 



We now replace the Lorentz group by general linear transformations, considering Lagrangians with an affinc 
symmetry. The fields now transform as 

&<P = ?<P,i+\e\{v<P) k i, (18) 

with £ l = a 1 + e l k x k , with general (constant) e i k , and where the action of the general linear group on the fields 
is given by 

(a(p)\ = scalar (19) 

{aAi)\ = 25? Ai vector (20) 

(crhi m )\ = 2(5fh im + 5%hu) tensor. (21) 

Assuming invariance of the action means that C transforms as scalar density and thus, SC = (££ 1 )^ = C^C + 
j. For the rest, the argument goes just as in the case of the Lorentz symmetry, with the only difference that 
we end up with 16 (instead of 6) independent equations. The result is 

T \ — ~2^ m k,m (23) 

This is now in full analogy to equation Q, in the sense that invariance under linear transformations (together 
with the global translations) dictates the form of the conserved current t\, which is again in the form of a 
divergence. 

In contrast to the Lorentz case, where the Minkowski metric 77^ had been introduced by hand, no metric is 
needed for the above arguments. However, in order to make contact with the results obtained for the Poincare 
group, let us define 

Then, we have obviously S mlk = L m ^ k \ and from the expressions (|16f) and J23J, we can derive the following 
expression for the symmetric Belinfante tensor 

rpik ^£7nik j^ikm _|_ j-^imk _|_ j^rriki _|_ j^kmi j-^kim^ (25) 

A few remarks are in order at this point. First, if T l k has the form of a relocalization term, i.e., if L m \ in (IL'MI) 
is antisymmetric in im, then it follows from (|25l) that T lk = identically. Second, it should be noted that 123(1 
and l|25(l are very strong relations. Indeed, for a scalar field, e.g., they lead (in view of l|19Jl ) immediately to 
r l k = 0. The reason, however, is also obvious: It is not possible to construct a theory that is invariant under 
general linear transformations only with scalar fields. For gauge fields (i.e., vector fields), we find 

r\ - -{^At). m . (26) 

We see that r l k i — is identically satisfied if dC/dAf m is antisymmetric in im. This is indeed the case in the 
physically relevant situations, where the derivatives of the gauge fields enter only via the Yang-Mills tensor F£. 
However, there could be exceptions, in particular concerning derivative couplings, e.g., of the tetrad field ef in 
gravitational theories without independent connection. This turns out not to be the case though, as we will 
now show in general. 
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4 General covariance 



From the explicit form (|23[1 . we see that only vector and tensor fields give explicit contributions to t\. Thus, 
the question is what kind of restriction do we get on the form of those contributions from general covariance of 
the theory? Well, this is not difficult to find out. Under general coordinate transformations x % — ► x % — £ l , the 
vector and tensor fields transform as 3 

SAt = CAi^+C^Ai, (27) 

Sgim = Cmm,i + C,l9im + C,m9ti, (28) 

while for the Lagrangian, we must have 5C = (^C)^. Similar as in the the case of gauge symmetry considered 
in the Introduction, it turns out that it is actually enough to consider transformations of second order in x\ 
i.e., = s l kl x x . In any case, for the physically important case, namely the vector fields, we find 

(i,m) 

or, in other words, 

L m l = -L im k . (29) 

The same result is found for the case of tensor fields, but this is not really of interest, since there is no generally 
covariant action for a tensor field containing only first derivatives. In any case, we will give a general proof (for 
theories containing up to second derivatives) in section 6. Equation l|29|) is the analogue of equation and 
simply means that the canonical stress-energy tensor l|23|) is identically conserved. In particular, the momentum 
vector 



P h = I Ado* = / r° t d 3 



.r 



can be written as a two dimensional surface integral 



P k = - j> L m \da mi =-- j L*%dV (30) 

Moreover, as mentioned before, a consequence of H29|l is that the Bclinfante tensor vanishes identically (a result 
that has been found in 0]). 

As a result, we have shown that in generally covariant theories, the canonical stress-energy tensor can be 
written in the form of a relocalization term, and thus, the momentum vector in form of a two dimensional 
surface integral (see also [5]), while the Belinfante tensor is identically zero. This holds for theories based on 
Lagrangians containing only first derivatives of the fields. 



5 Applications 

5.1 Special relativity 

The physically relevant theories that are generally covariant contain necessarily the gravitational field. Never- 
theless, in order to illustrate our formalism, let us consider the following Lagrangian which is independent of 
any metrical background 

3 Those relations are simply the Lie derivatives of the corresponding fields. 



G 



This Lagrangian is a total divergence, and therefore, the field equations are identically satisfied for any Aj. The 
stress-energy tensor, according to © is of the form 

i _ nliqp A A _ si p raqp a a 
' k ~ zt -^-q^p^l.k " k t Slrs^-pq, 

while from (JS3J, we find 

i _ njiqp A A 
I fc — z & -^q.p-^k^ j 

which is obviously identically conserved. Also, the difference of both expressions, e m9P Ff. m F pq — \8 k e lmpq Fi m F pq 
is easily shown to be zero. A similar calculation shows that the Belinfante tensor T lk is zero. 



5.2 Einstein-Cartan theory 

Einstein-Cartan theory is based on a tetrad field ef and a Lorentz connection T ab i . (Latin indices from the 
beginning of the alphabet a, b, c . . . refer to tangent space, with Minkowski metric rj a b — diag(l, — I , — f , — 1).) 
For completeness, we consider the Einstein-Cartan-Dirac-Maxwell Lagrangian 

£ = -l«R a \k ~ \F ih F ik + e [^ 7 m D m V - Vmh m iP) ~ mW>], (31) 

withe = dete?, R ab lk = r^-F^+r^F^-F^L^, F lk = A^-A^ and = d^- {T a \G ab ^-iqA^. 
From lO and (gO) , we find 

t\ = -(ee^rr ab fc ), m - {eF* m A h \ m . (32) 

Obviously r l k i — identically, and T' tk = 0, according to our general theorem. The expression is easily 
generalized to non-abelian gauge fields Af. 

In order to evaluate the momentum P k in l|3U|) . it is useful to make a few general considerations. First, it 
is clear that, in the presence of radiation, we will, in general, not find a finite expression when the surface of 
integration is extended to spatial infinity. On the other hand, if we assume that the leading order contributions 
have the behavior Ai — 0(-), and thus Fi k = 0(-i?), then it is clear that the second term in i.Y3 does not 
contribute to P^ 4 . The same holds for non-abelian gauge fields, and in particular for T a \ itself, if we include in 
C terms quadratic in the curvature R ab ik - The only exception to this is the tetrad field, which can be assumed to 
behave like ef = Sf + 0(~). Thus, quite generally, the only contributions that enter explicitely the momentum 
P k are those stemming from the Einstein-Cartan Lagrangian (first term in l]32[l). and eventually from additional 
terms in the Lagrangian, quadratic in the torsion tensor T? k = e£ i + r° bi e^ — e a ik — F a bk e b . Such terms, however, 
will in general also modify the Newtonian limit of the gravitational theory. 

Note that for the same reasons, relocalization terms are usually considered not to modify the momentum 
P k , and are therefore used as a tool to modify the stress-energy tensor according to our will, without changing 
the physically important quantity, which is P k . If this were generally true, however, then in generally covariant 
theories, the momentum would always be zero, since we can always write the stress-energy tensor in the form 
of a relocalization term. It is therefore very important that, as we have argued, the linear Einstein-Cartan term 
eR provides an exception to this. On the other hand, for the same reason, it is quite a questionable procedure, 
in the framework of those theories, to modify the stress-energy tensor by relocalization terms, and, e.g., define 
the Belinfante tensor. Those procedures were invented in special relativistic theory exactly because they are 
supposed not to modify the momentum, and not simply to overcome our dissatisfaction with an asymmetric 
tensor. After all, as we have shown in the introduction, the electric current density too can be written in the 
form of a relocalization term, i.e., j* = C lk k , with antisymmetric C %k . However, would anyone ever come up 

4 For the boundary conditions necessary to obtain a finite expression for the conserved charge, see also 0. 
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with the idea to add an additional relocalization term because he does not like the specific form of j l ? And 
even if he did, he would certainly take care that at least the charge itself is not modified by this procedure. 

Let us return to the expression (|32|) . Recall that the connection can be split into a torsionless part and 
the contortion tensor T"^ = + K a \, where, in Einstein-Cartan theory, the contortion tensor is directly 
expressed in terms of the spin of the spinor field. In the absence of spinor fields, we have K a \ = 0, and we can 
evaluate T ab i using ef k + T a bk e b — efT l ik , where T l ik are the Christoffel symbols. For the Schwarzschild metric 
ds 2 = (1- M/r)dt 2 - (1 - M/ry 1 dr 2 -r 2 (dd 2 +sin 2 M(p 2 ) and using a diagonal tetrad, we find P k = 2ttM<5£, 
i.e., the energy is just one half of the mass, Pq — m/2 (M — 2km, with 8nk = 1 in our units). The same 
result was found in the preprint version of [3j ■ Note that the same result emerges from the Reisner-Nordstroem 
solution (since only the leading order term of goo = — l/ffrr = 1 — M/r + q 2 /r 2 contributes to the surface 
integrals at spatial infinity). 

Since r l k is asymmetric, and T' lk — identically, we cannot formulate a conservation law for the angular 
momentum based on the stress-energy tensor. However, we can exploit the Lorentz gauge symmetry of the 
theory, i.e., the invariance of the Lagrangian under 

Se*=e a b e b , ST ab t = -D 4 e afc , 6i/> = - l -e ah a ab ^ : (33) 

where e ab = -e ba and D,e ab = e ab . l + T a cl e cb + T b cl e ac . Let us define the following current 5 



, - _ 1 BC 1 dC i 8C dc dC r . 

ab - 2 del k ^ 2 de b k &al 4 c^, fc ^ dT cb i k at 5T» C . fc 6i ' ' ' ' 4 1 



which is antisymmetric in ab and can be interpreted as angular momentum density current. Performing the 
same steps as in the introduction, we find from global symmetry (e ab constant) 

J ab k ,k - 0. (35) 

and from linear symmetry (e ab i k = 0) the explicit form 



Jab" = (g^b-h, (36) 



and finally, from local symmetry that the expression in parentheses is antisymmetric in ki, such that the current 
is identically conserved and the corresponding charge can be put into the form of a two dimensional surface 
integral. Note that for the particular case of Einstein-Cartan theory, the first two terms in l|34|l are automatically 
absent. 

Although we focused, in this section, mainly on Einstein-Cartan theory, it is obvious that the results are also 
valid for any candidate of Poincare gauge theory, i.e., for any theory with an independent Lorentz connection 
and with a Lagrangian at most quadratic in curvature and torsion. We refer to [S], where the Hamiltonian 
analysis of those theories has been carried out. The analysis is easily extended to metric affine theories (with 
an independent general linear connection, see [7]), provided one finds a way to couple the spinor fields to the 
connection, see, e.g., 

Finally, we should also mention that in some sense, the Lorentz gauge group is more closely connected to 
the coordinate transformation group than conventional internal gauge groups. In order to perform Wigner's 

5 For simplicity, we omit again the terms with <tjj. 



8 



classification of elementary particles, you need to analyze the behavior of the fields under both the diffeomor- 
phism and the Lorentz group. The link is provided by the fact that in the flat limit, ef = Sf, and the residual 
transformation freedom is then, apart from global coordinate translations, a global Lorentz gauge rotation and 
a simultaneous global Lorentz coordinate transformation (with the same parameters), such that ef = Sf re- 
mains unchanged. In other words, the Poincare transformation group of special relativity (on which the particle 
classification is based) emerges in the flat limit as a combination of both the Lorentz gauge group and the diffeo- 
morphism group. As a result, it is not unnatural to consider combinations of both groups right from the start. 
That is, instead of the conventional Lie derivatives Sip — <fi,iC + k( a( P) k i related to pure diffeomorphism 
invariance, one can consider modified Lie derivatives under which the theory is still invariant, as a result of the 
additional Lorentz gauge invariance. In this way, one can derive in quite a natural way alternative expressions 
for the conserved energy- momentum currents. The details of this procedure have been recently worked out in 



5.3 Tetrad gravity 

A more conservative way to incorporate spinor fields into general relativity is by simply replacing the metric 
tensor with the tetrad field, without introducing an independent connection. The gravitational field then couples 
to the spinor field via Bitp = d l ip - where T a \ is a function of the tetrad field and its derivatives that 

can be evaluated from e° fe + r a bfe ej = efT l ik . The free gravitational Lagrangian is taken in the form 

£ = -l[-\r m r ikl - l -r^r lk% + \r m lm T l \\ (37) 

where T a ik = eg f — e" fe . This Lagrangian coincides up to a surface term with the Hilbert-Einstein Lagrangian, but 
it has the additional feature that it does not contain second derivatives of the tetrad field, and it is generally 
covariant under coordinate transformations. Note also that C is invariant under (local) Lorentz rotations 
5ef = e a h e\. For this Lagrangian, we find 

L mi k = e[r fc ™ + r im k - r im k - Sir 1 "} - 5^r% (38) 

which is antisymmetric in im, and thus leads an identically conserved T l k and to T lk = for the Belinfante 
tensor. 

Note that, as a result of the derivative couplings, in the presence of spinor fields (with the same Dirac 
Lagrangian as in H31(l . where Dj^> is replaced by Diip) we find an additional term in the form 

T mi(D) jr rni „ im ( Jml /QO\ 

L k = 4 Ffc ~ a k + cr k\i ( 39 ) 
where c Qb m = is the so-called spin density of the spinor field. (This is not a conserved quantity, though.) 

Again, i l m ' l ^ D ^ leads to a surface term. (There can also be additional terms, from gauge vector fields, in the form 
of the second term in (|32|1 which, as we have argued before, do not contribute cxplicitely to the momentum.) 
The theory has a local Lorentz symmetry 

6e? = e a b el 6j> = - l -e ab a ab ^ 7 (40) 

which could eventually be used to define an angular momentum density for the system. However, from the 
above symmetry, we find relations formally identical to QH4JI. (|35|l and (|36|l . without the terms involving T ^. 
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In particular therefore, from (|3fcil> . we find that the current vanishes identically. Thus, in this theory, it is not 
possible to formulate a conservation law for angular momentum, neither from the stress-energy tensor, nor using 
the local Lorentz symmetry. Note however that you can use the vanishing of the expression l|34|) in order to 
simplify the stress-energy tensor obtained from (|38|) and (|39fl . by observing that the second term in i|34|l is again 
the spin density o a k . We then find L m ^ = —2a lkm (for the total expression, plus and therefore, 

from l|15|) . r[ ,fe l = a' lkm m , showing that in the absence of spinor fields, r lk is automatically symmetric. 



5.4 General relativity 

Finally, we consider classical general relativity. Since we have restricted our formalism to Lagrangians that do 
not contain second derivatives of the fields, we start with the Lagrangian 

(41) 

which is equivalent to the Hilbert-Einstein Lagrangian — h ^J—gR. Apart from gik , there might be gauge vectors 
A" and scalar fields ip, which do not contribute explicitely to Pk (and do not appear at all in the Belinfante 
tensor T lk ). We assume that there are no derivative couplings of gik to any other fields, as is indeed the case with 
gauge fields and minimally coupled scalars. (Spinors cannot be described within a purely metric formalism.) 
After some algebra, we derive from l|41ll the following expression 



j^mk 



Im ck kl cm . mq Ik kq h. 

9 ,i d i ~9 ,A +9 i9iq9 ~9 i9i q g 



+ [V=9 (-2g lm S k + g lk 5? + g km 5\)] , . (42) 



The expression in the first row is antisymmetric in mk and is thus not only identically conserved, but moreover, 
it does not contribute to the Belinfante tensor T lk , as we have established previously. The expression in the 
second row, however, is not antisymmetric in mk. This is because C is not covariant under general coordinate 
transformations. It is not hard to see that nevertheless, the divergence of this expression vanishes identically. 
The expression is in the form of what we will call in the next section a second order relocalization term. Also, 
since £ is covariant under linear transformations, the relation r l k — ~\L m \. m is still valid. Note, however, that 
in order to evaluate the total canonical stress-energy tensor, you have to add the contributions from the vector 
fields to L m \- (Although they will in general not contribute to Pk anyway.) 

Finally, since L mk is not antisymmetric in mk, the Belinfante tensor T lk does not identically vanish, and 
indeed, we find from 

Tik 1 f / f^ik Im im kl km li , iklm\~\ //10^ 
= 2 W~9 KV 9 -g V ~9 v +g V )\ t i im - ( 43 ) 

This is exactly the form derived in [3] by the direct application of the Belinfante relation to the canonical 
tensor © (and a subsequent use of the Einstein field equations) and the tensor T %k is known as Papapetrou 
tensor. Note that matter fields (vector and scalar) do not contribute to T lk , since they lead to antisymmetric 
(and zero) contributions in L™, and therefore T lk represents the total stress-energy tensor of the system. 

It is interesting to remark that T' lk , although apparently not written in the form of a relocalization term 
(see, however, section 6), is identically conserved and it is also not hard to show that the momentum can again 
be written as a two dimensional surface integral, although two steps have to be performed to achieve this. (In 
the integral over T ok from l|43|l . split first one of the indices I or m into time and space components, and then 
the other one. The term with both I = m = does not contribute, while the remaining terms lead to surface 
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integrals.) It turns out (see that the same Belinfante tensor emerges from the full (covariant) Hilbert- 
Einstein Lagrangian. (Not the same canonical tensor r l k , though.) This shows that the statement given at the 
end of section 4 (i.e., that in generally covariant theories, the Belinfante tensor vanishes identically) does not 
generalize in its full extend to generally covariant theories with second derivatives in the Lagrangian. (However, 
there is a generalization of that statement, in the sense that the stress-energy tensor r l k is still identically 
conserved and the momentum can be expressed in terms of two dimensional surface integrals, similar as is the 
case with T lk here. We will show this in section 6.) 

We have evaluated the momentum emerging from t\ and T lk and find P k = AttMS^ = m5° k for the 
Schwarzschild metric in both cases, which is the double value of that found in Einstein-Cartan theory. This is 
in accordance with the results found in 3 , where in addition, the tensor T lk has been shown to be equivalent, 
as far as the leading order contributions are concerned, to the corresponding Landau-Lifshitz tensor and the 
Weinberg tensor. (Note that only the leading order contributions contribute to the surface integrals at infinity 
anyway.) 

The Belinfante tensor (|43|l being symmetric, we can formulate a conservation law for the angular momentum 
density T lk x m — T lm x k , which would not be possible from t\ or from any other line of argumentation (as was 
the case in Einstein-Cartan theory). 

As a final remark, it is important to have in mind that the fact that the stress-energy tensor (both r l k and 
T lk ) is identically conserved can not be derived merely from the symmetry of C under affine transformations, 
but it emerges here rather as an additional, unexpected result. Otherwise stated, not every Lagrangian covariant 
only under affine transformations leads to an identically conserved current. The specific Lagrangian considered 
here, however, has the additional property that it is equal, up to a four divergence, to a generally covariant 
Lagrangian, albeit one of second order. As a result, under general coordinate transformations 5x l — £ l , the 
Lagrangian transforms as 5C — + i- e - ; h picks up an additional four divergence. Quite obviously, 

such a transformation behavior is sufficient for the field equations to be generally covariant. 

5.5 Discussion 

Summarizing our results, we found that in Einstein-Cartan theory, the canonical tensor is identically conserved, 
while the Belinfante tensor vanishes. Nevertheless, it is possible to establish a conservation law for angular 
momentum, based on the local Lorentz gauge invariance of the theory. Thus, in a sense, there is no need for 
a symmetrized stress-energy tensor. The separation of the angular momentum from the stress-energy tensor 
(which, after all, is the current corresponding to translational invariance) seems to be satisfying also from a 
logical point of view, if we recall that the spin structure of the spinor field was completely separated from the 
coordinate transformations (i.e., ip is treated as a scalar field), and was transferred to a tangent space with an 
inherent Lorentz symmetry. Thus, it should also be that same Lorentz symmetry that is responsible for the 
angular momentum conservation. (The argument is not absolute, though, since after all, the angular momentum 
is not entirely given in terms of intrinsic spin.) 

On the other hand, in general relativity, we have a non-vanishing symmetric Belinfante tensor, and we can 
therefore formulate an angular momentum conservation law based on this tensor. Nevertheless, we see two 
problems with this procedure. First, in the context of a generally covariant theory (or at least, covariant under 
the general linear group, in the first order theory), we arbitrarily pick out a specific subgroup, namely the 
Lorentz group, and modify the canonical stress-energy tensor, which was already in the form of a relocalization 
term, by a relocalization term (of the same order of magnitude) such that the resulting tensor allows for the 
formulation of the conservation laws corresponding to that subgroup, namely angular momentum conservation. 
This, in our opinion, is not quite in the spirit of general relativity, which is based on the full diffcomorphism 
invariance, and has no inherent preferred subgroups. (Again, the argument is not absolute, since nevertheless, 
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the Poincare group (or eventually the de Sitter group, in presence of a cosmological constant) emerges as 
groundstate symmetry of the theory, if the correct signature of the metric is assumed a priori.) Second, the 
metric framework of general relativity does not allow for the presence of spinor fields. Therefore, as soon as 
we deal with spinor fields, we have to go over to the tetrad formulation, which has therefore to be considered 
as more fundamental. However, as we have established, the tetrad formulation leads to a vanishing Belinfante 
tensor, and there is no obvious way to define a conserved angular momentum current. Therefore, the apparent 
success of the Belinfante tensor in general relativity seems to be only a coincidence that does not generalize to 
the more fundamental formulation of the theory. And most importantly, the absence of a conserved angular 
momentum in the tetrad formulation of general relativity provides an argument in favor of Einstein-Cartan 
theory with independent Lorentz connection. 

Our opinion is therefore that the use of the Belinfante symmetrization procedure should be confined to 
the purpose it was initially designed for, namely to special relativistic theories, where the Poincare symmetry 
is an inherent ingredient right from the start, and where the integrated momentum is not influenced by the 
relocalization terms as a result of the asymptotic behavior of the fields in conventional theories. Both of those 
requirements are violated in general relativity. 

Nevertheless, in the framework of general relativity, the tensor (|43H can certainly have its usefulness, if 
interpreted correctly and most importantly, if used consistently. It represents, after all, a conserved current 
and is, in this aspect, no different from the Weinberg or the Landau-Lifshitz tensor (see |3|). But it is also 
clear that any expression of the form (|43|l . with g lrn replaced with any other symmetric tensor, rf m replaced 
with any other symmetric tensor, and \J—g replaced with any other function, is identically conserved too, is 
symmetric too, and is also related to the canonical tensor by a relocalization term (of second order, see next 
section). Thus, modifying the canonical tensor by relocalization terms and requiring Poincare invariance does 
in no way fix the form of the stress-energy tensor. 



6 Second order field theory 

For the sake of completeness, we now extend our analysis to theories based on Lagrangians containing second 
derivatives of the fields (we refer to such theories as second order). The second Noether theorem tells us that 
again, invariance under a local symmetry leads to an identically conserved current. We will establish the explicit 
form of the current corresponding to general coordinate covariance. 

As can be expected, the analysis for second order theories contains one more step, namely global, linear, 
quadratic and finally general gauge or coordinate transformations. We start again with gauge theory as a 
warmup exercise, but since we do not know of any physically relevant second order gauge theory anyway, we 
confine ourselves this time to the abelian theory. 

Thus, we consider a theory invariant under 

5Ai=E ti , STp = -ieip. (44) 
Having in mind that the field equations (for a generic field ip) in second order theories are of the form 

d£ , d£ . . d£ . 

),*» + («, ),m,i=0, (45) 

dip d<p, m 0<P,m,l 

and requiring S£ — under ()44[). first with constant e, then with linear e, then with second order e (i.e., 
e = £ikX l x l with constant e^) and finally with general e(x), we arrive at the following four equations 

J*.i = (46) 
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dC dC dC 



= 
= 

where the current has been defined as 

J 1 = 



dC 



i,k,l 



dC 



dC 



dC 



iip 



PL 

dip.ij 



■ifp,k + ( 



dC 



(47) 
(48) 
(49) 

(50) 



Note that in 149(1 , symmetrization over the three indices is understood. The situation is essentially the same as 
in first order theories (see equations (0, 10, (@J). Global symmetry (|4"fi|) leads to current conservation, linear 
symmetry 1(4 7JI enables us to write the current in the form of a divergence, and finally, local symmetry (in this 
case, quadratic and cubic) shows that the current is identically conserved. More specifically, ((481) shows that 
the two first terms in ((47(1 are identically divergence free, while l(49|) shows the same for the last term in ((47(1 . 

Note that, if we take into account 1(4811 and 1(49(1 . then the first two terms in (14 7() are in the form of the 
previously encountered relocalization terms, C %k k with antisymmetric C lk . It is interesting that the last term, 
although identically conserved, is not of the same form, but rather in the form C tkl k ; with C %kl such that the 
totally symmetric part is zero, C^ %kV> — 0. It is not hard to show that such a term in the current J % leads again 
to a two dimensional surface term for the charge J J°d 3 a;, namely (n, A = 1, 2, 3) 



G 



ha 



,0,fc 



d?x = 



C , fc d <7 A 



a 



<0k0 



d J a„ 



where we use the fact that C 000 = 0. Therefore, we will refer to such terms as second order relocalization terms 6 . 
Note that the stress-energy tensor 1(43(1 consists of such a term. 

We now turn to coordinate transformations in second order theories (see \W\). Thus, we consider again 
transformations of the form 

8<P = ?<P,i + l?,kW\, ( 51 ) 

where the expressions for (cr(p) k i are given in l(19|) . I(2U(1 and 1(21|) . Then, we require invariance of the action, i.e., 
5C = (C^).i and derive the corresponding conservation laws. The manipulations are simple and we give only 
the results here. 

From invariance under global translations (£' = a 1 = const), we find the conservation law for the canonical 
stress-energy tensor 

r\ k = (52) 

with 

- S k C-(- \ lt p. 



dC 



Invariance under general linear transformations (£' = e' 1 k x k , with constant £ l k ), determines again the form of 
the stress-energy tensor 



-<P,i,l- 



(53) 



3 The quantities C k and C nkl are also referred to as superpotentials. 



dC 



.1 



dC 

d<p, m ,i 



(54) 



,m,Z 
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Invariance under second order transformations (£ l = £ikX l x k ) leads to 



dC 



1 dC 



dp>,k,: 



2 dp., 



dC 



dpj 



[(°Y>)*],J 



(55) 



{km) 



where the subscript (km) means that the expression in brackets has to be symmetrized in km (no differentiation) . 
Finally, the requirement of general covariance (£ l (x) arbitrary) leads to 



dC 

9<p, m , 



0, 



(56) 



(krnl) 



where the expression is totally symmetrized in kml (no differentiation). Quite obviously, (|55|l shows that the 
three first terms in l|54|l are in the form of a relocalization term C kl il , with C \ = —C \, while (|56[1 shows 

that the last term has the form of a second order relocalization term C klm ilm , with C^ klm \ — 0. Thus, r\ 
is identically conserved and Pk can be written as a two dimensional surface integral, in accordance with the 
second Noether theorem. 

From l|54|) . it is also clear why the Belinfante tensor cannot vanish in second order theories. That is because 
the Belinfante symmetrization procedure is based on the definition (|16|l . More precisely, the Belinfante procedure 
for second order theories (see |3]) consists in writing (|54|l in the form 



k _ jj_ t mk 

i 2 i,m> 



with 



L 



m k 



. d£ 

2^ <P, 



dC 



(apf 



2 J^-\{aip)\ 



Qp,k,m ' Otp.m 1 dtp.j^m 

and then defining (compare (|16|l ) the Belinfante tensor 



dC 



5" 



D J, 



(57) 
(58) 

(59) 



where S mkl — and L mkl — rf l L m f. This tensor is easily shown to be symmetric if the theory is Poincare 

invariant, i.e., if at least the antisymmetric part of (|54|l is satisfied. Moreover, if the theory is invariant under 
the affine group, i.e., if H54|) is satisfied, then T tk can be written in the form of equation (|25|) . i.e., 



T 



i k 



" ikm 



L 



/ink 



L 



mki 



L k 



knii 



'- kim~\ 



^j^(m,i)k j^{ik)rn _|_ j^(rnk)i^ 



(60) 



If the theory is generally covariant, then we find from l|55|) for the symmetric part of L mkl the following relation 



]_{mk)i 



dC 



dp, 



dC 

dip,j,k 



(apY 



(61) 



from which you can explicitely evaluate T lk . It does not identically vanish, but it is nevertheless identically 
conserved, as a result of 15611 . In the case of general relativity, with the covariant, second order Lagrangian 
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£ = kV~9-R*i it leads again to the expression lj4"3T) . as is easily verified, in accordance with the results of 3 7 . It 
is interesting to remark that the evaluation of T lk involves a lot less computation than that of r l k . 

We wish to point out, however, that you can always write the last term in l|54|l in the form of a first order 
relocahzation term too. The has been shown in |l(Jj for the general case. Specifically, in general relativity, this 
term is equal to 

-i [^=5 (-2g lm S k + g lk 6? + g km Sl)] ^ , 

which corresponds to the second line in (IS) (hence the identical Belinfante tensors in first and second order 
theory). An equivalent form of this expression is 

which is now of the form — ^L lk i ; , with L lk i antisymmetric in kl. Together with the three first terms in 1)54(1 . 

we can now write the stress-energy in the form t\ = —^L i h with antisymmetric L i.e., in the form of a 
first order relocahzation term (although still equal to <|57|) h There is no reason to prefer one or the other form, 
meaning that the generalization of the Belinfante procedure to second order theory is not really as unambiguous 
as it has been presented in 3 . Indeed, defining the Belinfante tensor using L in l|60() leads to a vanishing 
tensor. 

Finally, we note that (|61fl reduces in the case of first order theories to (|29H . which completes the proof of 
the results of section 4 without explicit reference to the vector or tensor nature of the fields. 

The above results do not mean that there is any fundamental difference between first and second order 
theories. In both kind of theories, the stress-energy tensor is given in the form of relocahzation terms and 
can thus obviously be annihilated by adding appropriate relocahzation terms, of first or second order. Also, 
in both cases, relocahzation terms modify the momentum vector in the case of the linear Hilbcrt-Einstein or 
Einstein-Cartan type theories. The only difference comes from the following: In first order theories, if you 
eliminate the antisymmetric part of T tk , making it thus suitable for the formulation of conservation laws for the 
full Poincare group, then you automatically annihilate the complete tensor, while in second order theories, you 
can eliminate the antisymmetric part and still have a non- vanishing tensor. However, as we have argued before, 
in generally covariant theories, it is rather arbitrary to pick out the Poincare group and modify the stress-energy 
tensor accordingly. There seem to be only two reasonable (albeit extreme) points of view: Either, we allow for 
modifications of r l k by relocahzation terms or we don't. If we do, then we should use them in order to construct 
the tensor that has the largest possible symmetry. This is quite obviously the tensor that vanishes identically, 
since it is trivially fully covariant. If we do not, then we should simply stick to the canonical tensor in its initial 
form. Everything else is, in the framework of generally covariant theories, an ad hoc procedure and can only be 
of limited usefulness. (See, however, the remarks at the end of section 5.) 

7 It is understood that the matter part of the Lagrangian is first order, and thus does not explicitely contribute to T lk . If this 
is not the case, you will have, apart from 1431 . additional terms in T tk , which, however, do not contribute to P k if the usual 
asymptotical behavior is assumed. A similar assumption has been made in [3], where the Belinfante procedure was applied to 
the canonical tensor of the gravitational field, and then brought into the form 1431 by the use of the gravitational field equations. 
The result is then the Belinfante tensor of gravity plus the Hilbcrt tensor of matter. However, the Hilbert tensor is not equal to 
the Belinfante tensor of the matter field if the matter Lagrangian contains second order derivatives (see next section). The only 
generally valid form for the total Belinfante tensor T lk is therefore I6UI with 1611 . 
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7 Matter stress-energy and the Hilbert tensor 



In this section, we split the Lagrangian into a gravitational and a matter part, C = Ca + C rn , and briefly examine 
the relations between the canonical stress-energy for the matter fields, the corresponding Belinfante tensor and 
the so-called Hilbert tensor. Since the relation between canonical and Hilbert tensor has been discussed in detail 
in our previous article and since the generalization of the discussion to the Belinfante tensor is rather simple, 
in view of the relations we have already obtained in the previous sections, we will confine ourselves to give a 
very short exposition of the issue. The subject has also been covered, e.g., in an d 

We consider purely metric theories (e.g., general relativity), and for (some kind of) completeness, we allow, 
for the moment, derivative couplings (first order) of the metric to the matter fields. Nevertheless, non-minimal 
couplings of the form y/— g ip 2 R, are still excluded, since they contain second derivatives of the metric coupling 
to ip. Moreover, we confine ourselves to first order theories, i.e., we assume that no second derivatives of the 
matter fields occur in the Lagrangian. 

The variation of the matter Lagrangian then is of the form 



SC m = 



dip 



-Sip - 



dip, 



dgir 



-Sg, 



lm 



dgirr, 



Sg 



lm,i i 



(62) 



where ip denotes collectively the matter fields (i.e., all the dynamical fields, except for gi m )- Note that the 
matter fields are supposed to satisfy the field equations, while for gi m , no field equation can be used, since the 
gravitational part Cq is not included in our Lagrangian. In other words, gi m is treated as background field. 
Under a coordinate transformation, the fields transform according to Q51[). and in particular, for the metric 
we have Sgi m = gi m ,k^ k + £ fc ;fffcm + £ k . m 9ik- The matter Lagrangian is assumed to be a scalar density, i.e., 

The Hilbert stress-energy tensor is defined by 



Sgir, 



= -2 



( dC r 



dg, 



! Ill 



Irn.i 



(63) 



In view of the 



where, in accordance with our previous conventions, we have defined T lm as tensor density, 
invariance of C m under coordinate transformations, T lm satisfies the covariant conservation la' 
Y™T lk = 0. The derivation of this law is found in any textbook on general relativity (or, see ^J)- Note however 
that it holds independently of the specific form of the gravitational Lagrangian. 

With this definition, and using T'™ = 0, we derive for £ l = a 1 the following conservation law 



X _ rr% 
k,i k,; 



dgin, 



0, 



(64) 



where r l fc is defined as in © (but with C m ), and T\ = T tm g m k- Next, for linear transformations, we get 



n+-L 



mi _|_ cy 



dg 



li.m 



-9kt 



3£ m 
dgim, 



■gim,k — 0, 



(65) 



where L™, is defined as in l|22|l . Requiring general covariance, we find the additional relation 



- (mi) 



j dC'm 
dgi(i,m) 



9ki = 0, 



(66) 
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showing that t\ — T\ + ® g f m . gim,k is identically conserved. It is not hard to see from the definition i|16|) that 

the Belinfante tensor is not equal to T jfc , nor does it reduce to T tk in the flat limit — rjik- However, in the 
case where derivative couplings are absent, l|6*5|) simplifies to r l k — T\ = ~\L ml k m , and, since reduces to 

l^ mi k = 0, the Belinfante tensor T lk is now easily shown to be directly related to T lfc , namely, we have 

T ik = T l m rj mk , (67) 

where we recall that T\ n = g m k?' lk , and thus, since the second index of T' lk is lowered with the Minkowski 
metric, we can equivalently write 

T\ = T\. (68) 

In summary, in first order theories without derivative couplings, the Belinfante tensor is indeed equal to the 
Hilbert tensor (when written in mixed form). Note that the Belinfante (matter) tensor T lk is, strictly speaking, 
not symmetric (since it is not equal to T lk in contravariant form). A symmetric tensor is obtained by raising 
the second index of T\ with g %k instead of rf k . This, however, is not really of importance, since T % k is not 
conserved anyway in a curved background (neither is T tk ), and in general, neither will the angular momentum 
be conserved. 

The equality of the Belinfante and the Hilbert tensors in first order theories in the framework of metrical 
theories of gravity can give rise to two different interpretations. The first one is that it provides a strong 
argument to the standpoint that the Belinfante procedure can indeed be used to derive (or rather define) the 
correct stress-energy tensor not only in special relativity, but also in metric theories of gravity, and in particular 
for the gravitational Lagrangian itself. A second, quite different point of view consists in interpreting the Hilbert 
tensor as a covariant generalization of the Belinfante tensor, a point of view equally strongly supported by the 
above equality. According to the second interpretation, the stress-energy tensor of the gravitational field itself is 
then given by the variation of Cq with respect to the metric, and quite obviously, the total stress-energy tensor 
would vanish in view of the gravitational field equations. This interpretation is supported by the fact that in 
generally covariant first order theories, the total Belinfante tensor is indeed identically zero, while in second 
order theories (like general relativity), one could modify the Belinfante symmetrization procedure, including 
second order relocalization terms, in a way that again, the total tensor would vanish. 

We have carried out a similar analysis of the relation between the canonical and the corresponding (general- 
ized Hilbert) tensor T\ = — in the framework of Einstein-Cartan theory and also in tetrad gravity in |llj . 
It was shown in particular that both tensors are not in general related by a relocalization term (because of the 
additional field V ab i in the first case, and because of the derivative couplings to spinor fields in the second case), 
and it is not hard to extend the analysis to include the Belinfante tensor. The result is that there is no equality 
(not even up to a relocalization term) in those theories between the Hilbert and the Belinfante tensors. For 
instance, in Einstein-Cartan theory, one can derive the relation T l k — T l a e k — o- ab l T ab k , with a ab l — 9£ m /9r ah i . 
Equality holds in the absence of spinor fields a ab l = or in the limit of vanishing gravity, ef — Sf,T ab k = 0. 
Similarly, in tetrad gravity, we have equality only in those limits. 

Moreover, the above result for metric theories does not generalize to matter Lagrangians containing second 
derivatives. This is actually quite obvious, since if it would, then the same would be true for the gravitational 
part of the Lagrangian, and then, the total Belinfante tensor (adopting the definition of [jj], i.e., using 1581) 
and 1)60(1 ) would be equal to the variation of the total Lagrangian with respect to gik, and thus zero, which is, 
however, not the case in second order theories, as we have shown in the previous section. The equality between 
Hilbert and Belinfante tensor can be achieved by modifying the Belinfante relation for second order theories, 
as we have done in the previous section, rewriting the second order relocalization term in the form of a first 
order term, such that the total Belinfante tensor vanishes again. However, in view of the physical irrelevance 
of matter Lagrangians with second derivatives, we do not pursue this idea further. 
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8 Generators and Hamiltonian constraints 



In this section, we briefly review the relation between the conserved charges, the generators of the corresponding 
symmetries and the (first class) Hamiltonian constraints. The issue is of fundamental importance for the 
quantization of field theories with local symmetries and has been pioneered by Bergmann and Dirac. We 
confine ourselves to a brief analysis of the important cases (internal gauge symmetry and general covariance) 
and refer the reader to the initial articles, in particular ^3], ^H], [T^], ^7] and |18|. 



8.1 Gauge theory 

In order to present our equations in a covariant form, we use the covariant formalism of |19j . In particular, we 
will use integrals over spacelike hypersurfaces a defined by <I>(x) = 0, with (timelike) normal vector = The 
normal vector is understood to be normalized, g tk riink = 1. Note that this normalization is only a convention 
used for greater convenience in intermediate steps. Our initial and final relations should not depend on this 
normalization, since we do not want to assert any special meaning to the metric guc, which is considered to be 
a field like any other field. 

Further, we define the canonical momenta by 

T (i) = |^ (69) 

where we put the index into parentheses to remind of the fact that only the part normal to the hypersurface, 
7r = 7T^n, corresponds to the physical momentum known from conventional Hamiltonian theory. 

We start again with non-abelian gauge theory. The conserved current density has been derived in the form 



Next, in view of equation Q we define 8 



and in view of equation , 



Fir 

Ja ~ dAf k %a ' [U) 



as well as the corresponding charges 

Q = f £ a Jida h Q = I e a r a d^, Q = I [e^^da,, (73) 



where e a (x) is an arbitrary parameter. 

As a result of equation we have the identically satisfied relation J l a + J l a = 0, and as a result of 
3 % a + Jq* = 0. It is not hard to recognize that the corresponding equation Q + Q = is equivalent to the 
secondary (first class) constraints of the theory. Consider, e.g., flat space and choose = <5° (that is, the 
hypersurface t = const). Then, the integrand of Q + Q reduces to e a [— i-K^a a tp + irgC^Aj. — n£ ], where 

8 Note that this is the opposite of the expression if GJ is taken into account. 
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fi = 1,2,3. In the last term, we have exploited the antisymmetry of %a resulting from Q. This is exactly 
the secondary constraint arising in Yang-Mills-Dirac theory. In the context of quantization in the Coulomb 
gauge, it is usually referred to as Gauss' law. More generally, the relation J l a + J* = is equivalent to the field 
equations for the Yang-Mills field, but the fact that it can be derived directly from the symmetry of the theory 
indicates the presence of constraints. The field equations for %j) for instance cannot be derived in this way. 

In a similar way, the relation = is related to the primary constraints. Consider again the case n, = 5®. 
Then, we have 7r^ = 71^ , and therefore 7r° = 0, which is the well know primary constraint of Yang-Mills theory. 
Alternatively, with the help of equation (|111|) of Appendix A, we can express the primary constraints simply as 
P = 0. 

Next, we analyze the relations of Q, Q and Q to the generators of the gauge transformation. For this, we 
assume the following commutation relations between fields and canonical momenta 

fc>(s),¥>(y)]=0, [*®(*U (m) (}/)]=0, y{x)^{y)]=5i{x-y), (74) 

where x 7 y are assumed to be separated by a spacelike distance. For details on our formalism, we refer to 
Appendix A. In the case of spinor fields, the above relations are replaced by anticommutation relations, e.g., 
{i/j m (x), 7rjy (y)} — 8^5 l (x — y), where M,N denote the spinor components. It is not of our concern here 
where the (anti)commutation relations come from, be it from a classical Poisson bracket (see or be it from 
postulating them in a quantum theory. They are simply assumed to be valid. 

We are now ready to evaluate the commutation relations between the fields and the charges defined in 
For the manipulations involved, a few helpful rules are provided in Appendix A. Here, we give only the 
results. Note that in order to derive those results, two dimensional surface terms have been omitted during the 
calculations. This is permitted if we assume that the parameter e a has an appropriate asymptotical behavior, 
or if we simply assume that it vanishes outside of a certain region. 

We find the following relations 

ie a <r a ip, (75) 
0, (76) 

0. (77) 

In order not to overload the notation, we have omitted the arguments of the fields. It is understood that all 
the fields (including e a ) are taken at a point y, lying on the hypersuface with respect to which the charges have 
been defined in (J73J, except in the second terms of (|?6*j) and iffTjl. where e a is taken at the point x, and the 
integration is performed over do~i(x). As pointed out in Appendix A, there is no obvious way to simplify those 
integrals without specifying a hypersurface, but they do not appear in the combinations of the generators we 
are interested in. If we choose the hypersurface t — const for /j, — 1,2, 3, we find the more conventional form 9 

[Q,Al\ = -e% [Q,A?] = -e% 

9 Note that the same results can be directly derived by the use of the more familiar, but not covariant generators 

Q = j e a J°d 3 x, Q = J e a J°d 3 x Q = J [e a j°°}, d 3 x, 

instead of 1731 . where we identify directly the physical momenta with n = tt' ' , and take the commutation relations to be (at equal 
times) [ip(x),ir(y)] = 5^(x — y). 



[Q,Af] = s a cjA], [Q,ip] = 

[Q,Af] = -e^+e%n m n l+ J e^SKx-y)]^, [Q,1>] = 
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as well as [Q, A%] = [Q, A@] = 0. 

Those results are very interesting. First, we see that the charge Q does not (as is occasionally stated) 
generate gauge transformations on the fields. Rather, it generates homogeneous transformations on both Af 
and ip, in the adjoint (or vector) representation on Af (a rotation in isospin space) and in the fundamental 
representation on ip. 

The inhomogeneous part is generated by Q + Q, namely we have 

[Q + Q,Af] = -e^. (78) 

In order to find the generator of the full gauge transformation, we need all three parts, i.e., 

[Q + Q + Q, Af] = -e? - e^ a A], [Q + Q + Q,tp] =ie a a a ip, (79) 

which is exactly the transformation Q we started from. Both 178( 1 and (|79|) are independent of the hypersurface. 
Another combination of interest is the generator Q + Q, because it is identically zero (related to the secondary 
constraint). If we choose the hypersurface t — const, we find for the spatial components [Q + Q,Af] — 
— e@ — e" „ a Al, The fact that an operator that vanishes weakly (i.e., as a result of the constraints) generates 
gauge transformations on the propagating fields (consider, e.g., the Coulomb gauge), is interpreted in quantum 
theory as the expression for the fact that physical states have to be singlets under the gauge group. That is, free 
bosons Af cannot be part of the physical particle spectrum (see, e.g., [201 ) - This should be somehow alerting, 
considering the fact that we will now go over to spacetime symmetries. 

8.2 General covariance 

The analysis is quite similar to the previous case. In view of the relations 10, (|23(1 and (|29|) . we define 10 

r\=^^ :k -SlC (80) 

r\ = Af\ m = -^k ( ' 1 W"], m , (81) 

and 

r m i=\L mi k = \^ m \^)\. (82) 
Next, we introduce the corresponding charges 11 

-k ~i j _ £> / r^fc 



P = Je k r\da t , P = j£«f\d<j t , P = / [e*f m 4], m do-i. (83) 

The primary constraints P — arise again from the fact that T m l is antisymmetric in mi, while the secondary 
constraints are expressed by P + P = 0. The explicit form will depend on the nature of the field, i.e., on the 
form of (aip) 1 k . Note also that the relation r l k + f % k — is equivalent to a field equation. In general relativity, 
this would be the Einstein equation. In other theories, it could be a combination of field equations. 

In order to evaluate the commutation relations with <p, it is important to recall that in (|80|) , the Lagrangian 
is expressed in terms of the field and its derivatives. It is not understood that field derivatives are to be replaced 

10 Note that again, the second expression 1811 corresponds to the opposite of 1231 . 

11 Again, there is the more familiar choice P = J e k r° k d 3 x, P = J e k f° k d^x and P = J [e k T 00 k ] ,od 3 £, involving only the compo- 
nents 7n ' of the momenta, which satisfy [ip(x),Tr(y)] = <5 3 (x — y). 
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by momenta. This cannot be done in a unique way anyway, since (|80|) is not a Legendre transformation. In 
other words, T l k has nothing to do with the field Hamiltonian, which is indeed a Legendre transformation (see 
|T§1 for more details and on the definition of the Hamiltonian in the manifestly covariant formalism) . 
We find the following commutation relations 

[P,tp\ = -ip, k e k , (84) 

[P, d = -\e k m {°vTk + jIMV],^ + \ J [fax ~ vliWk^&r™ ( 85 ) 

[PM = -\[(*v) k m e m hn k n* - \ J [fax - y)M<T<p) m k e k da m . (86) 

In particular, for the sum P + P, we have 

[P + P^]=-\e k m {^n- (87) 

The situation is in complete analogy to the previous case. The charge P (i.e., the canonical field momentum 
in the strict sense) generates what is usually referred to as spacetime translations, in the sense that it tells 
us the evolution of if from one point to another, (p(x) — tp(x') = —ip y k£ k - On the other hand, the operator 
P + P generates passive coordinate transformations as they are usually considered in general relativity, i.e., 
<p(x)—<p'(x') — —^e k m {af) m k . And finally, the operator P+P+P generates the active coordinate transformations 
we started from, i.e., the Lie derivatives of the field 

[P + P + P, ^] = -<p, h s k - \e k m {ovT k = ~ <P&), (88) 

with x n = x l + e l . 

A few remarks are in order at this point. The relation H84J) is in the form one expects for a momentum 
operator. It is the operator one wants to have in a quantum theory. And it is based on the canonical stress-energy 
tensor. Now, since in the framework of general relativity, T l k +f % k — is equivalent to Einstein's equation, i.e., to 
SC/Sgik = 0, the tensor corresponding to the gauge generator P+P is obviously given by —yJ~gG lk +T' lk , where 
G lk is the Einstein tensor and T lk the Hilbert tensor (density) for the matter fields. In particular, concentrating 
on the matter part, we see that T lk , which is also equal to the Bclinfantc (matter) tensor T lk , is related to a gauge 
transformation of the form (for the hypersurface t = const) Sip = —s k tp : k — | £fe m( <7 ¥')"fc + ^[( ai P) k £k ]fi- This 
makes it rather hard to interpret the quantity J T° k d 3 x as field momentum. Quite obviously, this interpretation 
should be reserved to the corresponding expression with the canonical tensor. 

In short, the Belinfante is not directly related to the generators of coordinate transformations and as such, 
its integral over space should not be directly interpreted as energy and momentum. 

A second remark concerns the same generator, P + P, which is zero (secondary constraint), and generates 
again gauge transformations on the propagating fields, quite similar as in the case of internal gauge theories. 
For instance, for the spatial components of a tensor (e.g., the metric) we find (for the hypersurface t = const) 
[P + P, g^v] = —g iaJ ,k^ t — £ ,u,9ku — £ k v 9^k- For transformations restricted to the hypersurface, we have e k = 
(0, e^), and thus, we find that [P + P, g^ v \ is equal to the Lie derivative of g^ v in the three dimensional subspace. 
One is thus tempted to conclude that physical states have to be singlets under coordinate transformations, i.e., 
scalars. This would not only exclude gravitons, but also photons and vector fields in general. While for the 
latter, we can argue that they are considered usually on a given background (with a reduced symmetry), in the 
case of the gravitational theory, the solution to this problem is more profound and can be found in the fact that 
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general relativity (or similar theories, e.g., Einstein-Cartan) is actually a theory with a spontaneous symmetry 
breaking, in the sense that the vacuum is given by go? = (or eventually a de Sitter metric) and thus allows 
only for Poincare (or de Sitter) transformations. This invalidates the previous arguments, and gravitons can 
now be interpreted as fields propagating on this background. (Note that similarly, massive gauge bosons in the 
Salam- Weinberg model are also part of the physical spectrum.) 

Finally, we caution that we have assumed in our analysis that the theory is free of second class constraints. 
Second class constraints arise in particular in the Dirac theory, but also in Einstein-Cartan theory (see |Hj) and 
are not related to gauge symmetries. Their presence leads to modifications of the canonical (anti)commutation 
relations and to modifications of the above results, in particular of the relation 184|) . We refer to f° r details 
on the Dirac theory. As to Einstein-Cartan theory, one has to take into account both the diffeomorphism 
covariance and the local Lorentz invariance H33[) . This leads to the primary constraints (identifying directly 7r 
with 7r(°) for simplicity) Tp ah = tt® = 0, where TT l ab and 7r* are the momenta canonically conjugated to and 
ef, respectively. The remaining primary constraints, namely = and n£ b + |(e°e^ — e^e°) = 0, are not 
related to any kind of symmetry and are second class, i.e., the commutators [tt", ir^ + | (e°e^ — e^e°)] do not 
all vanish, as can be seen by expressing ee^ in terms of ef . 



8.3 General relativity 

The previous analysis is valid for generally covariant first order theories, and as such, is directly applicable to, 
e.g., Einstein-Cartan theory or Poincare gauge theories in general, once we have dealt consistently with the 
second class constraints. On the other hand, general relativity is either second order, or the Lagrangian is not 
generally covariant and requires therefore special care. We will not deal with the second order theory here, but 
instead, try to fix the first order theory in order to make it suitable for the above analysis. 
For simplicity, we consider pure gravity, such that gik is the only field, and we have 

jmk = d*L {agin)ki = A J^- gu = AaPWgu. (89) 



dgin 1 dg, 



ki. 



The explicit expression is given in (|42l) . In particular, as a result of the non-covariance of the Lagrangian (|41|1 . 
£mk - 1S n0 ^ an tisymmetric in mk. The relation l|89[) is actually quite interesting, since it tells us that there are 
no first order generally covariant theories (i.e., based on a covariant Lagrangian) for symmetric tensor fields at 
all. Namely, since n kl ( m ) is symmetric in kl, it cannot be antisymmetric in km, except if 7r fei ( m ) is zero (that 
would be the case for a Lagrangian consisting only of a cosmological constant). It is not hard to derive the 
same conclusion for scalar densities (so-called dilaton fields) and for symmetric tensor densities. 

On the other hand, in order to evaluate the stress-energy tensor, we are not really interested in L m \ , but 
rather in L m \ m . Therefore, as we have already indicated in section 6, we can modify the second line in l|42(l and 
write instead \\f Z7 g{— 2g lm &\ + 2g lk 8 r j a )},i- We denote by L m \ the expression (f^ with the second line replaced 
in that way. This is now antisymmetric in km. Then, we have 

~ L rak _ L rnk = [^/=g^kgm _ ^™^)] ; . (90) 

Obviously, L m \ m = L m \ m . We can thus write t\ = —^L m f m , with L m \ antisymmetric in mk. In view of l]89p. 
we introduce modified momenta ft 11 ^ by requiring 

L^=4jf fe| W fflj , (91) 

which leads to 

~ki(m) ^ki(m) , f t kl\ mi ( t km\ li /nn\ 

* - t y ' + W-99 ),i9 - W-99 ),l9 ■ (92) 
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Note that 7f fci ( m ) is no longer symmetric in ki. Next, we postulate canonical commutation relations between 
fields and modified momenta in the form (for spacelike separations) 

[g ik (x),7r ln ^(y)} = \{5 l M + 5 l k m 6™(x-y). (93) 

Only the symmetric part of 7r'"( m ) does not commute with gi k and contains thus the physical momenta. 
The stress-energy tensor r l k — (d£/dgi rn .i)gi m ,k — 5 k C can now be written in the form 

r\ = n lm ^ gim>k - 5\C = n lm ^g lm , k - [{V~gg p[ ), P g vm - (V~gg u )^g pm ]gi m , k - SIC. (94) 

We can now proceed as in the previous cases, i.e., we define the generators by 

P = Je k T\da u (95) 

as well as 



P = / e k f \da t (96) 



with t\ as above and with 



Finally, for the third generator (see (JS3J), we have 



k = - 2 L k,m = ~ 2 F 'gki\,m. (97) 



P = Jler m k l }, m da l = 2 J [e k n^], n d<j u (98) 

where f™ = ^L m k . The primary constraints are now expressed by the fact that f™ is antisymmetric, and thus 
P = 0, while the secondary constraints are given by P + P = 0. Thus, in a sense, with the modification of the 
canonical momenta, we have restored the primary constraints which have been lost because we started with a 
non-covariant Lagrangian. 

The commutation relations of those generators with the metric lead consistently to the expressions (I84L 
H85|l and 186JI . Explicitely, if we choose again the hypersurface t = const, we have 

[P,9ik] = -g lk , m e m , (99) 
[P,g lk ] = ~e m l g mk -e m k g mi + (e m g mk )^ + (e m g mi ) fi Sl (100) 
[P,9ik] = -(9im£ m ),o5 k -(9kme m ), Q 6l (101) 

From this, we can see explicitely that action of the constraint P + P = on the spatial components, namely 

[P + P, g^} = -g^,x£ X - e x ^g KV - £ A ,. 9km , (102) 

where we have assumed that e° = 0, i.e., the transformation takes place on the hypersurface in question. As 
pointed out earlier, this is again the Lie derivative of the three dimensional metric. In Appendix B, we analyze 
the question whether <7 M „ can indeed be assumed to represent the propagating part of gi k in an appropriate 
gauge. 

In order to check whether the introduction of the new momenta fully consistent, one has to introduce 

a Hamiltonian H = J T\n k dai (see ^H])> express it in terms of the momentum eliminate 
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the velocities gi r 



and analyze the commutation relations of H with the momenta. This is a non-trivial 



task, in particular one will have to deal with ordering problems, see, e.g., It is however not unusual in 

the context of first order general relativity to carry out modifications by hand, in order to achieve consistency 
between the constraints of the theory. For instance, Dirac [22] chose to modify the Lagrangian by a surface 
term in a way that the momenta conjugated to go^ vanish weakly, which then represents a primary constraint. 

An alternative way, and probably a more elegant one, is to start from the covariant, second order Lagrangian. 
Let us briefly sketch how this could work. One defines the canonical momenta 



dips 0(p tk ,i 
and writes the stress-energy tensor (|53f) in the form 



P 



m(i) 



dC 



(103) 



(104) 



where rp m = ip m plays the role of the canonical variable conjugate to p m w. The identically conserved form of 
the stress-energy tensor from (|54|) takes the simple form 



1 r 
2 



(105) 



where {p^j) k i — 2Sj(p,i + [(cy) i\,j> i- e -> ^ ac ^ s correctly on rpj in accordance with its total tensor structure, 
taking account of the additional vector index. Similarly, we express equations (|55|) and H5(jf) in terms of the 
momenta and finally construct four charge operators, the first two expressing the constraints stemming from the 
equality between l| 104(1 and (|105fl . and the other two expressing the vanishing of the expressions (|55|l and (|56|l . 
For instance, in a metric theory, it can easily be seen from l|56ll that one set of (primary) constraints is related 
to the fact that no second time derivatives of the the components got are contained in the Lagrangian. We 
thus see that this apparent coincidence in general relativity is actually a necessary feature of generally covariant 
second order theories. 

The only non vanishing commutators in second order field theory are given by 



b0r),7r«(y)] =SUx-y), [<Mar),/«(y)] = 8 k Mx - y), 



(106) 



where x and y are assumed to be separated by a spacelike separation. 



9 Conclusions 

We have analyzed Noether's theorem localizing step by step the transformation group from global to local. The 
intermediate steps of linear (and eventually quadratic) transformations turned out to be of importance, not 
only in general relativity, where one usually prefers the use of a first order Lagrangian which is covariant only 
under affine transformations, but also in theories based on a fully covariant Lagrangian. Indeed, the restrictions 
on the conserved current derived from each step of locality are directly related to the first class Hamiltonian 
constraints of the theory. The primary, secondary, etc., constraints related to the invariance of the theory under 
a certain symmetry group, as well as the corresponding generators, can be directly read off from the form of 
the Noether current as it arises in each step. 

Concerning the stress-energy tensor, we have shown that it is identically conserved whenever the theory is 
generally covariant and the momentum can thus be written in terms of two dimensional surface integrals. In first 
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order theories, the symmetric Belinfante tensor was shown to be identically zero. In second order theories, this 
is not necessarily the case, but the generalization of the Belinfante formula to those theories does not seem to 
be unique. Special attention has been paid to general relativity, which, in the first order approach, shares many 
features with theories based on a fully covariant Lagrangian, but at some points, modifications are necessary in 
order to restore the properties that have been lost as a result of the use of a non-covariant Lagrangian. 



Appendix A: Covariant Hamiltonian formalism 

We recall the main features of the manifestly covariant formalism for field theory used in |19| . 

Let x = x 1 = (x° , x 1 , x 2 , x 3 ) be spacetime coordinates such that a hypersurface element can be written as 



d<Ji(x) 



( dx 1 (a)dx 2 {a)dx i {a) \ 
dx (a)dx 2 (a)dx 3 (a) 
dxVjdxVjdxV) 

\ dzVjdzVjdzV) / 



where dx'(a) means that dx l is restricted to some hypersurface a defined by $(2;) 
a; = const, we have da; = and do~i(x) — <5°d 3 x.) 
In the same coordinate system, we define 



5\x-y) 



f 5 x 0y0 


Six 1 


-y v )5{x 2 


~y 2 )6(x 3 


~y z ) 
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- y°)S(x 2 


-y 2 )S(x 3 
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-y :i ) 
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-y 2 ) 


) 
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0. (E.g., for the hypersurface 



(108) 



The transformation behavior for S l (x — y) under a coordinate change is found from the known transformation behavior 
of d<Ji(x) {y/—g d(7i is a vector) by requiring 8 l (x — y)do~i to transform as scalar under general coordinate transformations. 
Thus, 8 l (x — y) transforms as vector density. 

Next, consider a spacelike hypersurface a defined by = 0, with the (timelike) normal vector m = <!>,;. For 

convenience, $(1) can be chosen such that n 2 = riitikg = L Then, we have 



f{x)5\x - y)dai(x) = f(y) 



(109) 



where the integration is carried out over the hypersurface a containing the point y. For the specific hypersurface 
t = to = const, we find 12 , e.g., 



f(x)S t (x-y)dffi{x) = / f(x°,x)8 (3) (x~y) 8 x o y0 d 3 x = f(x°,y) 8 x o y o, 

J a 



(110) 



where x° is to be taken on the hypersurface in question, i.e., x° — to. Thus, if y° = to (i.e., if y lies on the hypersurface 
t — to), the result is simply f(y), while else, we find zero. Thus, 8 1 (x — y) can be seen as covariant generalization of the 
three dimensional delta function. 

Two useful relations are the following (see |19p 



f,ida k 



and 



/.fcdo";. 

/ f(x)n l do k (x) = / f(x)n k da t (x). 

J a J cr 



(ill) 

(112) 



12 We use both x° (and y° etc.) as well as t for the time coordinate, while to always refers to a constant. 
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Note that in general, Uj = rii(x), but we will omit the argument whenever there is no danger of confusion. In order for 
the above relations to hold, an appropriate asymptotical behavior of / has to be assumed, such that two dimensional 
integrals over the boundary of a can be omitted. 
In particular, we have 

/ f(x)S\x - y)A(n = I f(x)6 l (x~y)n k n k dai 

f(x)S\x~y)n i n k da k (113) 



which is equal to f(y) if y lies on the hypersurface. (We take the convention that all quantities whose arguments are not 
written explicitely are to be taken at the point x.) Let us introduce the following definitions 

da = nMfTj (114) 

8<r{x-y) = S'(x-y)m (115) 

8l(x — y) = 5 a (x — y)n = S m (x — y)n m n . (116) 

We can thus write 

f(x)d\x ~ y)d*i(x) = f f{x)5l{x - y)da l (x) = f f(x)So(x - y)da = f(y) (117) 

Jo Jo 

where for the last relation, it is assumed that y lies on the hypersurface. Moreover, we have ff{x — y)rii — 5 a {x — y)rn. 
Nevertheless, one should not confuse S l (x — y) which is given explicitely by (110811 . with 5%{x — y), which is defined with 
respect to a specific hypersurface. In particular, for t = to, we have S] r {x — y) = 5 x o y o 5^(x — y)g° l /g 00 - (The factor 
involving the metric components stems from the normalization n 2 = 1, which for t = to (and thus m = (no, 0,0,0) ) 
leads to no = g 00 and n l — g°'no.) In particular, in flat spacetime, we see that S^^x — y) has only one non- vanishing 
component, in contrast to (11081) . 

In order to derive the commutation relations with the charges we encounter in section 8, the following relation is 
useful 

/ f(x)[6t(x - y)],idffi(x) = -/,,(») + fAvWiy)^ + f f{x)%{x - y)],idat{x), (118) 

Jo Jo 

which is easily derived with the help of mill and 11121 . The last term in 11181 cannot be simplified without specifying 
a hypersurface. (Note that [S l (x — y)\,i = 0, but [5 a (x — y)\,% 7^ 0.) For the hypersurface t — const, we find by partial 
integration -/, M (g° M /3° )5", and thus 

/ f(x)\5i(x - y)],td<Ti(x) = -/,, + f.oSl (119) 

J o 

With the help of those relations, it is an easy task to evaluate the commutators 1751 - 177( 1 and 1841 - 187( 1. 



Appendix B: Propagating fields in linearized general relativity 

We consider a perturbation h ik on a flat background gik = r\ih + hit ■ The field equations of general relativity, \J~gGik = 
%k to first order in hik have the form 

h-Dtpik + i>\,k,i + 4> l k,id - Vik4> lm i, m ] = %k, (120) 

where ipik = hik — \h, and h — rf k hik- In particular, to this order, we have the conservation law T ,fc ^ = 0. Further, we 
have the gauge freedom of the linearized theory, 5hik ~ £i.k + £k,i, where is assumed to be of the same order as hik. 
All indices are raised and lowered with n^. 
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Our intention is to look for a Coulomb type gauge choice, where the non-propagating fields can be eliminated by 
solving the field equations (Gauss' law in electromagnetism) and the remaining fields satisfy a conventional wave equation. 
Note that, with the conventional gauge choice ip j = 0, the field equations take the form 7^ = — ^Otpik, which is easy 
to solve, but it does not tell us anything on the number of propagating fields (there are still 6 independent fields). 

A more convenient gauge choice can be read off the field equations after a 3+1 split, 

Too = itAVoo-^L) (121) 
To M = ^(Ai/>o M + Vv.o,!/ + ^"o iM ,k + V^o.^o) (122) 

= + ^o M ,o,* +^,0,^ + i> X u,\,» + i>\,\, v ) - ^v^(i' 00 ,o,o + i> XK \,K + 2^%,*), (123) 

where A = — 9 M 9 M and □ = 8q — A. We choose the following gauge conditions 

=0, ^ Om m =0. (124) 

The first condition on the three four- vectors ^ m is analogous to the Coulomb gauge in Maxwell theory, while the second 
one can be seen as a Lorentz gauge on the four- vector ip 0m . To show that it is indeed possible to impose such conditions 
on ip' k , consider the variations under gauge transformations of the quantities in question. We find 



Therefore, we perform first a transformation with £° satisfying D^ = —ip 0m m (and with £ M = 0), and then a transfor- 
mation with £ M satisfying A£ M = —ip^% (and with £° = 0). Both transformations exists generally, and since the second 
transformation does not modify the result of the first, the argument is complete. Note that it is not possible to require 
V> M "„ = together with = instead of 
The field equations take the simplified form 

%o = ^A^oo (126) 
= iA^ (127) 

Tflv = --£M>y.v + -(^On,0,v +lfov,0,ii - V^^° X ,0,x)- (128) 

The residual gauge freedom is now restricted to 

□£° = 0, + A£ M = 0. (129) 

Similar as in electrodynamics, we can eliminate the non-propagating modes by solving (I126t and 11271 . We have 

iMM) = --J^—^dx, (130) 

M,x) = _J_ /"^l^idV. (131) 

2tv J \x — x'\ 

It is an easy task to verify that this solution satisfies indeed the constraint ^} 0m m = 0. They are Coulomb type solutions, 
and for slowly moving matter distributions, the important contribution is given by l|130|l . corresponding to the Newtonian 
potential. 

To some extend, we have obtained a system quite similar to electrodynamics in the Coulomb gauge. As was to be 
expected, to the charge density correspond the four momentum density components %o, and the corresponding fields 
ipio can be eliminated by a Gauss type law. Also, just like A M M = 0, we have a remaining, symmetric transverse tensor 
satisfying tp^^ = 0, containing the propagating modes. Further, the non-propagating modes appear as an additional 
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source to the tensor ip^ in 11281 . This too is in straight analogy to the Maxwell equation OA^ = — j M + A°'^ in the 
Coulomb gauge. 

To simplify equation 112811 . we observe that at large distances from the matter distribution (far zone), we have 
approximately \x — x'\ ss \x\, such that the four equations (1130ft and 11311 can be approximated by 



ipoi{t,x) = ~ jiy / T 0i (t,x')d 3 x'. (132) 

Z7T \X\ J 

On the other hand, we have in the linear theory the relation T %k k — which leads upon integration to the conservation 
law {A/At) J T l0 d 3 x = 0. In particular, from l|1321 . we thus have ipoi.o = 0, and therefore equation 11281 reduces to 

□VW = -2T aP . (133) 

We also retain that from ipoi,o = 0, together with our gauge conditions l|1241 . we find 

if>",i = and ip " 3 ^ = for \x\ » £, (134) 

where I characterizes the dimensions of the matter distribution. Therefore, the conventional Lorentz type gauge condition 
holds again in this limit. Together with 11331 . we see that our field equations, as far as 7/> M „ is concerned, are completely 
equivalent to those of the conventional approach based on the gauge ip lk k . Nevertheless, we are one step ahead, because 
we did not simply forget the components ip^o, iftoo, but we have eliminated them properly via a Gauss type law. 

The above approach is probably the closest one can get to the Coulomb analogy of Maxwell's theory. There is only 
one disturbing point: The field satisfying the transversality condition i/> M ^ M = contains still three independent 
degrees of freedom, which is one more than those needed for the description of the massless spin 2 field. It is not hard 
to see that for plane wave solutions, the residual gauge freedom 11291 1 can be used to make traceless, thus excluding 
the existence of an additional scalar field. It should therefore be possible to eliminate a fifth component of i/) lfc either by 
an additional Gauss type law, e.g., the trace of tp lh or of ip^ v and relate it to a fifth charge (e.g., the trace of T lk or of 
T MI/ ), or simply to impose an additional condition of ip^ v to reduce the degrees of freedom to two. However, with the 
gauge freedom 11291 . this does not seem to be possible. Neither would the appearance of a fifth charge be physically 
satisfying. We therefore suggest that there exists a better gauge choice, such that upon eliminating four components by 
a Gauss type law and imposing four conditions on the remaining components, the gauge is completely fixed, and the 
degrees of freedom are reduced to two. We do not know whether such a gauge has been proposed in literature. 

In any case, we see that the propagating modes (in a perturbative approach) can be assumed to be contained in 
the spatial components of tpik, which is (up to a constant) the linear approximation of /_ g ffjfc, the tensor with inverse 

determinant of g ik . It is now an easy matter to evaluate the action of the generator P + P on this field, using dHD and 
(HJJ). We find 

[P + P, —j=g>j.v\ = -(^^£m^),ae A - {—^=g,>x)e\ v - {—=g u \)s x tll + (—^g^e^x, (135) 



= 0. This is the Lie derivative of the pseudotensor (1/^/— g)g^v, confirming our 
assumption that the secondary constraint P+P — generates coordinate transformations (on the hypersurface t = const) 
on the propagating components of the field, in complete analogy to the case of non-abelian gauge theory. 

Finally, we note that the gauge choice 11241 is the linearized version of {\f—gg 0l, ),i = and (V - gg^ v ),n ~ 0- 
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